Entanglement measures and the quantum to classical mapping 
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A quantum model can be mapped to a classical model in one higher dimension. Here we introduce 
a finite-temperature correlation measure based on a reduced density matrix obtained by cutting 
the classical system along the imaginary time (inverse temperature) axis. We show that the von- 
Neumann entropy S cn t of p^ shares many properties with the mutual information, yet is based on 
a simpler geometry and is thus easier to calculate. For one-dimensional quantum systems in the 
thermodynamic limit we proof that S en t is non-extensive for all temperatures T. For the integrable 
transverse Ising and XX Z models we demonstrate that the entanglement spectra of p& in the limit 
T — > are described by free-fermion Hamiltonians and reduce to those of the regular reduced density 
matrix pA — obtained by a spatial instead of an imaginary-time cut — up to degeneracies. 
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Entanglement is usually thought of as a quantum me- 
chanical entity yet it is well known that the properties of 
quantum models can be computed from a classical model 
in one dimension higher [IH3]- Such an approach is par- 
ticularly useful if one is interested in how entanglement 
builds up in a thermal state while lowering the tempera- 
ture or during unitary time evolution. 

Entanglement measures are, generally speaking, maps 
from the space of density matrices into the positive real 
numbers. Various entanglement measures have been in- 
troduced in the last decades and a set of axioms has been 
put forward which a good measure should possess [U[S]. 
For bipartite pure states the von-Neumann entropy of 
entanglement fulfills these criteria and provides a bridge 
into statistical mechanics and condensed matter physics. 
In particular, the reduced density matrix, whose eigenval- 
ues determine the entanglement entropy, is at the heart 
of the density matrix renormalization group (DMRG) [6J. 
This numerical method allows one to optimally approxi- 
mate pure states of many-body systems by matrix prod- 
uct states [7] and is most successful in one dimension. 

In experiments we are, however, usually dealing with 
mixed states about whose entanglement properties much 
less is known. Entanglement measures commonly used 
for pure states of multi-particle systems such as the en- 
tanglement entropy fail because they become extensive 
and thus no longer fulfill a boundary law. Furthermore, 
they cannot distinguish between classical and quantum 
correlations [4]. Entanglement measures which do dis- 
tinguish between these different types of correlations in- 
volve extremizations over all possible decompositions of 
the density matrix and explicit results have only been ob- 
tained for few-particle density matrices [SHIP]. Putting 
these fundamental difficulties aside, it is still useful to de- 
fine correlation measures for thermal ensembles which go 
beyond the one- and two-point correlation functions tra- 
ditionally studied in statistical mechanics and condensed 
matter physics. They might help to reveal, in particu- 



lar, phase transitions with complex or topological order 
parameters. 

In this letter we want to investigate how correlations, 
generated during imaginary time evolution, can be quan- 
tified. After a quantum to classical mapping we introduce 
as entanglement measure the von-Neumann entropy of a 
reduced density matrix obtained by a partial trace in 
the imaginary time direction. This measure also cannot 
distinguish between quantum and classical correlations 
but we will show that it shares many properties with the 
mutual information which has recently attracted consid- 
erable interest as a correlation witness for many-body 
systems |11H14| . The measure introduced in this letter 
is also a natural choice from the perspective of numeri- 
cal matrix product state algorithms: While the reduced 
density matrix obtained by a spatial trace is at the heart 
of the DMRG at T = 0, the reduced density matrix 
considered here is used in transfer-matrix DMRG algo- 
rithms to study the finite-temperature properties of one- 
dimensional systems in the thermodynamic limit [15-19]. 

The entanglement entropy for a bipartite system, S = 
A U B, is defined as 

Sent (A) = - Tr p A In p A (1) 

where pa — Tr# p, Tr pa — 1, is a reduced density ma- 
trix obtained from the density matrix p of the system 
S by spatially tracing out part B. If the system is in a 
pure state, then it is easy to show by a Schmidt decom- 
position that S en t = S on t (A) = S Qn t(B). It follows that 
S cn t cannot be extensive but rather has to scale with the 
surface between regions A and B [20]. In critical one- 
dimensional systems this boundary law is known to be 
weakly violated by logarithmic corrections |21l |2"2] . Here 
we want to study the mixed state described by the canon- 
ical density matrix p c = exp(~f3H)/Z where j3 is the in- 
verse temperature, H the Hamiltonian, and Z = Tr p c 
the partition function. Calculating the entanglement en- 
tropy for p c one finds that S cn t(A) S en t(B) in 
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general and that S en t becomes an extensive quantity ap- 
proaching the regular thermal von-Neumann entropy for 
P — > |21[ [23] . This can be easily understood as fol- 
lows: If all correlation lengths are much smaller than 
the extent of the considered subsystem then the rest of 
the system just acts as an additional bath. A way to 
partly correct this is to consider the mutual information 
given by I A>B = S cnt (A) + S ent (B) - S th (A U B) where 
Sth(^4 U B) = — Trpclnpc is the von-Neumann entropy 
of the whole system. The thermal contribution is then 
explicitly subtracted ensuring that Ia,b{P 0) —> as 
required. However, the mutual information is still not a 
proper entanglement measure because classical as well as 
quantum correlations contribute. Furthermore, the eval- 
uation of the mutual information even when using Renyi 
instead of von-Neumann entropies is quite involved due 
to the complicated geometry required to obtain general- 
ized partition functions [TT1 [12 HH HU HI] • An obvious 
question is if related finite-temperature correlation mea- 
sures can be defined which are easier to use in analytical 
and numerical calculations. 

Coupled qubits To motivate the idea of a correlation 
measure based on a quantum-to-classical mapping we 
start with two coupled qubits with Hamiltonian H = 
S\ ■ S2 where S is a spin-1/2 operator. Let us con- 
sider first the imaginary time evolution starting from the 
separable state \U) =_(\U)+ |+t»/2 + (|t4->- l+t»/2. 
Time evolving with e~P H this state becomes \^fp) = 

s£=^(rU}+ + CJ ^¥ M (\U)- lit}) and gets 

projected onto the maximally entangled singlet ground 
state for (3 — > 00. A picture showing how the entan- 
gled state is formed out of the separable state can be 
obtained by the quantum-to-classical mapping. We dis- 
cretize time into steps 5(3 writing p c = t m jZ where 
t = cxp(— 5(3 H) and M5j3 = j3. Rewriting the Hamilto- 
nian as H = iPj 2 — jldi 2 we find 



by defining a transfer matrix operator, a method of- 
ten used in statistical mechanics. One way to achieve 
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where Idi^ is the identity and Px,2 the operator per- 
muting the spins at sites 1,2 during the time step 5(3. 
Entanglement is thus being generated during imaginary 
time evolution by the braiding of the worldlines of the 
two spins. This is shown pictorially for one possible 
configuration of r-matrices in Fig. [TJa). For the sepa- 
rable initial state the correlations created can be 
measured by S cnt , Eq. Q, using p = |* ( 3)(* (3 |/(* (3 |* /3 ) 
as density matrix and tracing out one of the spins (see 
Fig.Qe)). When using p c as the density matrix , on the 
other hand, one finds S en % — In 2 independent of temper- 
ature. Thus S en t fails as a correlation measure and the 
whole temperature dependence of the mutual informa- 
tion, see Fig.[lje), stems from the thermal von-Neumann 
entropy with X A>B ~ 3/32T 2 for T -> 00. 

Here we want to pursue a different perspective on the 
classical representation of the qubits shown in Fig. |ljb) 
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FIG. 1: (Color online) (a) Local Boltzmann weig ht r, Eq. {2}. 
(b) Braiding of worldlines during imaginary time evolution 
discretized in steps 8(3 with green circles denoting tensor con- 
tractions, (c) Density matrix p c (l, 2|l', 2'), the lines connect 
contracted indices, (d) Density matrix p(l, 2, 3, 4|l', 2', 3', 4') 
projecting onto the ground state in auxiliary space, (e) En- 
tanglement measures as function of temperature. 

this is to perform an alternating 45° clockwise and 
anti-clockwise rotation of the r-plaquettes starting from 
Fig.JTJb). This leads to the completely equivalent graph- 
ical representation shown in Fig. [l|c) with plaquettes 
t r,l — Tr,iJ where Tr^ is the right/left shift opera- 
tor, respectively. Tracing either over the pair of open 
indices 1, 1' or 2,2' gives the reduced density matrix p A 
while tracing over both pairs gives the partition function 
Z. By moving the right 5 bonds (lines in Fig. [TJ C )) to 
the left and tracing over 1, 1' and 2, 2' we obtain the den- 
sity matrix p acting in auxiliary space along the imag- 
inary time axis, p, shown in Fig. [ljd), consists of two 
columns with the left column containing only right and 
left shift operators, T^l- Contracting pairwise the open 
indices 1, • • ■ ,4 on the left side with 1', •■■ ,4' on the 
right side yields again Z. We now investigate the re- 
duced density matrix p A — Tr^p obtained by taking a 
partial trace in auxiliary space. p A is then a 4 x 4 ma- 
trix whose eigenvalues can be easily calculated. Note 
that regions A, B always have to contain an even num- 
ber of Tfl^L-plaquettes so that the number of right and 
left shift operators is the same. A discretization using 
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more than four plaquettes would, in the simple case con- 
sidered, just add additional zero eigenvalues. The entan- 
glement measure S en t(A) = — Trp^lnp^ = S en t(B) is 
thus well-defined with lim^^o 5 ont = 2 In 2 and S cnt ~ 
3(1 + 6 In 2 + 2 InT) /64T 2 for T -> oo, see Fig. ge). 

Next, we want to generalize these considerations to 
one-dimensional quantum systems. A Hamiltonian with 
short-range interactions can always be written as H = 
^2jhjj + x, possibly in an enlarged unit cell. By us- 
ing a Trotter-Suzuki decomposition, the system can 
be mapped onto a two-dimensional classical system in 
much the same way we have mapped the qubits onto 
a one-dimensional classical system. A pictorial repre- 
sentation is shown in Fig. [2^a) where each plaquette 
is given by tr l — Tr,l exp(—S/3hjj^i). Note that 
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FIG. 2: (Color online) (a) Two-dimensional lattice obtained 
after the quantum-to-classical mapping. The shaded part is 
the QTM T . Performing the trace over A, B and the partial 
spatial trace over region B yields pA while tracing over A, B 
and taking the partial trace over region B yields p^. (b) 
Geometry to calculate Tr p\. The right part consists of two 
small cylinders while the left part is one big cylinder. 



metric representation is also possible but would require 
a wider column transfer matrix T ■ The reduced density 
matrix p^ — Tr§ p is then obtained by tracing over half 
of the indices along the imaginary time axis. Crucially, it 
is again easy to show using a Schmidt decomposition and 
choosing biorthonormal sets of right and left basis vectors 
that S en t = Sent(A) — S en t(B). p being a projector thus 
guarantees that S en t is non-extensive for all temperatures 
T. It is thus quite natural to replace the projector p onto 
the ground state, considered for zero temperature, with 
the projector in auxiliary space, p, for an infinite one- 
dimensional system at finite temperatures. 

Because pa is not a projector for T > 0, S en t(T) is 
difficult to calculate. Instead, the replica trick is often 
used [ST] which leads to the definition of Renyi entropies 
S n (A) = ln[Tx(//|)]/(l— n). A quantum-to-classical map- 
ping yields the geometry for Tr(p A ) shown in Fig. [2^b). 
Separate QTM's for the right and left part of the system 
can then be defined, however, an evaluation of S2 (A) re- 
quires one to calculate overlaps between the eigenstates 
of the QTM's [2] and is thus much harder than obtain- 
ing S^t. In the following, we study in more detail the 
scaling properties of S en t(T) and the spectra of p^ for a 
symmetric cut A = B = f3/2. 

Transverse Ising model We start with the transverse 
Ising model 

H = E (^K+i + ( 3 ) 

3 

where er are the Pauli matrices. This model shows a 
second order phase transition at A = 1 [3 . To study 
SentiT) and the entanglement spectra of we use a 
transfer matrix renormalization group (TMRG) algo- 
rithm |15ffTTl [15] . The entanglement spectra in the low- 



[hjj+x,hk,k+i] 7^ in general so that the mapping in- 
duces an error which is of order (5/3) 2 for the partition 
function. 

The canonical density matrix p c for periodic bound- 
ary conditions is shown pictorially in Fig. |2ja) with the 
l.h.s. and r.h.s. indices traced over. Tracing instead 
over the upper and lower indices and leaving the indices 
on the l.h.s. and r.h.s. open we obtain the density ma- 
trix p = T L jZ where T is a quantum transfer matrix 
(QTM) acting in auxiliary space. The partition function 
is given by Z = Tr^^ p c = Tr^S ?• For any T > the 
thermodynamic limit can now be performed exactly with 
limL^p = |*h><*l|- Here \% R ) and (V L \ are the left 
and right eigenvectors belonging to the largest eigenvalue 
of T with (^l\^r) = 1- This property can be easily 
understood physically because all correlation lengths — 
determined by the logarithm of the ratio of the leading to 
subleading eigenvalues of T — stay finite, p thus becomes 
a projector onto the ground state in auxiliary space. Note 
that in this representation p is non-symmetric. A sym- 
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FIG. 3: (Color online) TMRG results are denoted by symbols. 
Main: S cn t for various T and exact result for T = (solid 
line). Eigenvalues £j of for (a) A = 0.8 and (c) A = 1.2 
at T = 0.02. The free fermion levels are denoted by dashed 
lines, (b) S en t(T) for A = 1 and low-temperature fit S cn t = 
-0.1671nT + 0.402. 
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temperature limit are shown exemplarily in Fig. |3ja,b) 
and are equally spaced. The degeneracies of the dom- 
inant eigenvalues are given by 2, 4, 6, 12, • • • for A < 1 
and 1, 2, 1, 2,4, 4, 5,6, • • • for A > 1. As for the re- 
duced density matrix pa obtained by a spatial cut |25| 
these spectra can be explained by using corner transfer 



matrices (CTM's). We find p A = c 



-Hoi 
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where Hctm — ^2 e j n j is a f ree fermion Hamilton op- 
erator. For A > 1 all single particle levels are twofold 
degenerate, rij = 0,1,2, and given by €j = (2j + 
l)nK(y/l - l/\ 2 )/K{l/\) where K{x) is the complete 
elliptic integral of the first kind. For A < 1 we have 
e, = jnK(VT^X 7 )/K(X) where eo is non-degenerate 
and all other levels twofold degenerate. Knowing the 
spectrum in the zero temperature limit it is easy to cal- 
culate S cn t(T —> 0) |21| . The result is shown as solid line 
in Fig. [3] Right at the critical point the regular entan- 
glement entropy for an interval of length £ in an infinite 
chain shows a logarithmic divergence with system size, 
Sent = | lnf + Ci, with central charge c = 1/2 and a 
non-universal constant G\ |2T]. From this result it fol- 
lows immediately by a conformal mapping that 



S ent (T) = -Hv/T)+d 



(4) 



with v being the velocity of the elementary excitations. 
This result is universal for critical systems. For the trans- 
verse Ising model it is in excellent agreement with a two- 
parameter fit of the numerical data, see Fig.^b). 

XXZ model and boundary locality As a further ex- 
ample, we calculate 5 en t an d Pa f° r the XXZ model 



H = J 
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Here S is a spin- 1/2 operator and A parametrizes the 
exchange anisotropy. The model is critical for — 1 < A < 
1 with central charge c = 1. As shown in Fig. [4] for three 
different values of A, Eq. Q does indeed describe S en t for 
T« J. The spectrum of p A in the non-critical case A > 
1 can be constructed from a free-fermion Hamiltonian 
[2"5] and we find again that the same is also true for pa 
at low temperatures. As an example, the spectrum for 
A = 10 is shown in[4jb). The degeneracies of the free 
Fermion levels are the same as for the transverse Ising 
model with A < 1. however, now e 7 - = 2jcosh _1 A [25]. 
•^tot = 5Zj commutes with the Hamiltonian and is 
thus a good quantum number. For the QTM T it follows 
that S§. L = J2k(~^) k i^ii) r,l is a good quantum number 
where (S^)u t L acts a t site k in auxiliary space on the 
r.h.s. respectively l.h.s. of T. Matrix elements of T can 
thus only be non-zero if S z = — 5£ and eigenvalues 
of pa can be classified according to S z , see Fig.gb). For 
Pa it has been shown that its spectrum as a function of 
S* ot can be constructed for A ^> 1 by local perturbation 
theory in the spin exchange operator [26 . The same local 
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FIG. 4: (Color online) (a) One-parameter fits (lines) for 
T/J < 0.1 of TMRG data (symbols) according to Q with 
c = 1, v as known from the Bethe ansatz, and Ci as indi- 
cated, (b) Spectrum of p A for A = 10 and T/J = 0.02. The 
numbers give the multiplicities of each level (total in brack- 
ets), the dashed lines mark the free fermion levels. 



perturbation theory can also be performed for the QTM 
T. Keeping in mind that we have periodic boundary 
conditions along the imaginary time direction so that the 
partial trace introduces two boundaries, the degeneracies 
of each level for pa as a function of S z are exactly the 
same as for pa with two boundaries classified by Stof 
This demonstrates explicitly that a boundary law also 
holds for Pa obtained by a cut along the imaginary time 
direction. 

Ladders and two- dimensional models The quantum- 
to-classical mapping of a finite two-dimensional system 
of extent Li x L 2 leads to a three-dimensional system 
with dimensions j3 x Li x L 2 - In order to calculate the 
mutual information Ia,b following a spatial cut one can 
again use the replica trick. In particular, the Renyi en- 
tropy S*2 has been obtained from the three-dimensional 
generalization of the geometry shown in Fig. [2Fb) by 
quantum Monte Carlo (QMC) simulations fTTl 1121 12"4"]. 
However, the non-trivial geometry makes these calcu- 
lations rather complicated and 

Sent remains inaccessi- 
ble. On the contrary, the numerical calculation of S^nt is 
straightforward also in two dimensions. For L\ and L 2 
both finite, one can use the mutual information X A g — 
S en t(A) + S en t(B) + Trplnp to detect finite tempera- 
ture phase transitions. Furthermore, one can perform 
the thermodynamic limit in one of the spatial dimensions 
with \im Ll ^ 00 p(L 1 ,L 2 ,f3) = r) l\ becoming again 
a projector so that limT^oo S cnt = and a subtraction 
of a 'thermal part' is no longer required. 

Discussion We have analyzed the entanglement en- 
tropy S'cnt of a reduced density matrix obtained after 
a quantum-to-classical mapping and a partial trace in 
imaginary time direction. Our results show that trans- 
fer matrix DMRG algorithms are as efficient in simu- 
lating one-dimensional quantum systems in the thermo- 
dynamic limit at finite temperatures as regular DMRG 
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algorithms are for finite systems at T = 0. For two- and 
three-dimensional systems S cn t can be used to investi- 
gate finite temperature phase transitions. For QMC, in 
particular, this might be a viable and — due to the sim- 
pler geometry — easier approach than calculating Renyi 
entropies. 
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